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Buttiker-Thomas-Pretre (BTP) [Z. Phys B 94, 133 (1994)] proposed that the concepts behind 
the Larmor precession time tell us that it is possible to define exactly the local density of states 
(LDOS) in terms of the scattering matrix. However, we take into account evanescent modes and 
show that for an impurity in a quantum wire, this is in principle not exactly true. We also prove 
that the Wigner delay time gives correct superluminal times at the Fano resonances, in spite of the 
fact that the stationary phase approximation is not valid there. 

PACS numbers: 



I. INTRODUCTION 

Scattering phase shift in a scattering problem carry 
a lot of physical informations and are as important as 
scattering cross section or scattered intensity. However, 
phase shifts were always difficult to measure directly, un- 
til very recently 0, Q • Elastic scattering plays a promi- 
nent role in mesoscopic systems, wherein the inelastic 
processes are quenched by reducing the temperature and 
the sample size 0. Transport processes and thermody- 
namic properties of such samples can be formulated in 
terms of elastic scattering 0, Resonances in these 
systems, due to elastic processes, will be mainly Fano 
resonances and Breit- Wigner resonance is a special case 
of Fano resonance ^ . Impurities in such systems act 
like point scatterers and delta function like potentials are 
studied in this regard |^ Q . Such a negative delta func- 
tion potential can create Fano resonance. 



II. SCATTERING PHASE SHIFT AT FANO 
RESONANCE 

In this section we shall review some early works that 
lead us to study the present problem. Readers who are 
already familiar in this area can skip this section. 

Phase shifts for Fano resonance was first studied with 
respect to parity effect. Electronic states in a one- 
dimensional (ID) ring with an Aharonov-Bohm flux 
piercing the ring exhibit the parity effect according to 
which if the magnetization of the ring with N electrons 
is diamagnetic (or paramagnetic) then the magnetization 
with A'^ -|- 1 electrons is paramagnetic (or diamagnetic). 
Leggett conjectured that this is true even in presence of 
interaction between electrons as well as defects or disor- 
der in the ring |^ . Now suppose another quantum system 
S is coupled to the ID ring R such that the states of S can 
leak into R and become flux dependent. Ref. J^j shows 
that the states of the combined system (S-l-R) does not 
always have the property of reversing the magnetization 
with the addition of a single electron. It was found that 
the electrons in the ring undergo usual phase changes as- 
sociated with their quantum mechanical motion. These 



phase changes are i) Aharonov-Bohm phase, ii) statisti- 
cal phase, which means electrons being fermion, acquire 
a phase change of tt when they cross each other and iii) 
relative phase change due to the wave like property of the 
electrons, that depend on their wave vector or their ki- 
netic energy. It was also shown that apart from these 
phase changes, there are discontinuous phase changes by 
TT, at the zeroes of the Fano resonances (say at energy 
Ef)) that will be there when the ring (with the system 
S attached to it) is severed at a point and two leads are 
attached to the two broken ends 10]. If Fermi energy of 
N body state is below £'o, and that of the N + 1 body 
state is above Eq, then the magnetization of the N body 
state and the A^-|- 1 body state is the same, and otherwise 
opposite. Thus Leggett's arguments can be generalized 
to systems that are of the form S+R, with the conclusion 
that this new discontinuous phase change is a new phase 
different from the ones mentioned in i), ii) or iii) [lOj ]. 
When S and R are not ID systems but becomes quasi- 
ID (QID), even then this provides a general theory to 
understand coupled systems ITl|. 

In an experiment JJ, a quantum dot was coupled to 
a ring (S-l-R is now a dot-fring) and the conductance 
oscillations of the system with an Aharonov-Bohm flux 
was measured. Yeyati and Buttiker to interpret 

the conductance oscillations in terms of the flux depen- 
dence (or magnetization) of the electronic levels of the 
combined system (dot + ring) using the Friedel sum rule 
(FSR). The phase change at the resonances of the dot 
could be roughly understood but the phase change be- 
tween the resonances of the dot could not be understood. 
Other works tried to assign this in between resonance 
phase change to processes like charge decagture (the sys- 
tem throws away unit charge) [l^. Ref. ll and pre- 
dicted that the resonances are actually Fano resonances 
and the phase change between the resonances can be ex- 
plained by the new phase at the transmission zeroes and 
it will be an abrupt drop. It was confirmed in a later 
experiment that this phase drop occur over an energy 
scale, that is much smaller than any energy scale present 
in the system Finite width of the leads and evanes- 
cent modes has to be considered to explain the abrupt 
phase changes occurring between each consecutive reso- 



2 



nances. Ref. considers there is a pole between the 
zeroes and that is where a charge is captured and the res- 
onance phase changes smoothly by tt in agreement with 
FSR. But the discontinuous phase drops being a new kind 
of phase, does not have anything to do with charge cap- 
ture or decapture. Later on this was proved to be true 
whenever we have time reversal symmetry |l6l Il7j . Also 
the fact that simple ID calculations will not explain the 
phase drop between each two consecutive resonances was 
shown [l^l- Rigorous experimental verification has also 
emerged in favor of Fano resonances in quantum dots 

In fact unitarity is also req uired to produce such dis- 
continuous phase changes It was subsequently seen 
that when an unitary channel continuously evolves into 
a non-unitarity channel, then along with it, the trans- 
mission zeroes evolve into minima (difference with Breit- 
Wigner resonance is that the minima has a complex zero, 
that is for a complex value of incident energy, transmis- 
sion amplitude is zero), and discontinuous phase drops 
evolve into continuous phase drops |2ll |. As an exam- 
ple, one can consider a two channel quantum wire with 
a delta potential (see Fig. 1) at the middle of the quan- 
tum wire. Since the two channels of the quantum wire 
are opposite parity states like the ground state and the 
first excited state of a quantum well, the two channels 
are decoupled. If an electron is incident from the left in 
the first channel then amplitude of transition to the sec- 
ond channel i.e., ti2 (state 1 to state 4) or ri2 (state 1 
to state 2) is zero, making the first channel preserve uni- 
tarity. Whenever has a zero, the phase of drops 
(or decreases) discontinuously by tt. Now if the delta po- 
tential is shifted slightly from the center of the quantum 
wire, then parity is no longer a good quantum state and 
ti2 etc are no longer zero. So by continuously displacing 
the delta potential from the center, one can destroy the 
individual unitarity of each channel and thus make the 
discontinuous phase drops continuously evolve into fast 
continuous phase drops. These phase drops can be seen 
in Fig. 2 at an energy Ea « 85. Even dephasing can re- 
sult in the loss of unitarity [2^ 1 and can explain the small 
width of the phase drops observed in the experiment. For 
the two propagating channel case, 
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Here 6f stands for Friedel phase. For the delta poten- 
tial at the center of the wire, the discontinuous drops 
in arg{tii) or arg{t22) occur when |iii|^=0 and |t22|^=0 
and so it is obvious from Eq. 1 that there are no drops 
va Of. In that case, obviously the drops in arg{tii) or 
^^^5(^22) do not have anything to do with Of and hence 
charge decapture. But for the 5 potential off center, since 



1^11 P or 1^22 P are not zero when there are drops in their 
arguments, it is not obvious that these drops do not have 
anything to do with charge decapture. However only af- 
ter all the simplification, one finds that [2]| 
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That means Of only depends on the phase of reflection 
amplitudes and hence is completely independent of the 
phase drops in the transmission amplitudes. So once 
again the phase drops do not imply a drop vciOf or charge 
decapture 21]. 

One also knows from earlier known results that FSR in 
mesoscopic systems is valid only in semi-classical regimes. 
And so Fano resonance being a purely quantum interfer- 
ence phenomenon, one would expect larger violations of 
FSR at the Fano resonances. Instead what Refs. 
shows is that FSR is exact at the energy corresponding 
to the Fano resonance, but there are large violations away 
from this energy. Refs. also explained that the cor- 
rectness of FSR at the Fano resonance is due to the fact 
that there is a quasi bound state here, and hence the self 
energy due to the leads become minimum and hence its 
energy derivative becomes 0. This is obviously true for 
any potential that supports Fano resonances. The exact- 
ness of the FSR at the Fano resonance {Ea « 85) and the 
deviation away from it is also shown in Fig. 2 (compare 
the dotted and dashed curves) for the two propagating 
channel case. The results are the same for any number 
of propagating channels. The regime around Ea w 85 
is a purely quantum regime like the regime Ea < 40. 
However, in the regime around Ea « 85, which is the 
Fano resonance regime, FSR holds good unlike the regime 
Ea < 40, which is unexpected. 40 < Ea < 80 is the 
semiclassical regime and violation of FSR here is also 
unexpected. 



III. THE PROBLEM 

Apart from the DOS we know that time scales associ- 
ated with a particle crossing a quantum mechanical po- 
tential can also be determined from the scattering phase 
shifts. For example, in the the stationary phase approx- 
imation (i.e., phase shifts do not strongly depend on en- 
ergy) -^arg{ta/3) gives the Wigner delay time (WDT) 
for the particle to be transmitted from state a to state 
(3. A negative slope in arg(taa), like that in the solid and 
dash-dotted curve at Ea ~ 85 in Fig. 2, means super- 
luminality, i.e., the particle can travel faster than light 
across the potential, according to the WDT. However, 
once again, since Fano resonance is a quantum interfer- 
ence effect, dispersion will be very strong and stationary 
phase approximation will not hold good at the Fano res- 
onances, and one cannot be sure if these negative slopes 
actually correspond to super-luminality. To be sure of 
superluminality, established theories say that one should 
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see if we get negative delay times from the Buttiker- 
Thomas-Pretre (BTP) formalism. They proposed that 
the Larmor precession time can be determined exactly 
from the scattering matrix and give the correct local de- 
lay times in all regimes |^] and delay time can be de- 
termined by integrating the local delay times. It gives 



v(a,E,r, (3) 
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where v^a, E,r, (3) is proportional to the time spent by 
the particle (i.e., delay time) going from state a to state 
/3, while encountering the potential at r. If we sum it 
up for all a and f3 then it should give the LDOS at r, 
exactly 

Originally, it was derived by considering 
the effect of a small magnetic field on the outgoing spin 
wave function. But more generally, for any potential, to 
obtain the LDOS at r, we have to create a S function 
potential like local perturbation at r and see the change 
in the scattering matrix S of the entire system. Thus the 
delta function potential, apart from representing point 
defects, is also very ideal to study BTP formula and delay 
time. 

Given the fact that the scattering phase shift at the 
Fano resonance often violate established theories, is the 
BTP formula correct in QID? And secondly, is there su- 
perluminality at Fano resonance? WDT suggests that 
there could be limitless superluminality but can we trust 
WDT at the Fano resonance? 



IV. THE SCATTERING SOLUTION 

The scattering _2roblem is defined in Fig. 1. For this 
system in Fig. 1 
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where for the transmissions, a ^ (3 and 
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For a quantum wire with hard wall confinement, Pmn — 

{2m,-f/n^)sin[^{y, H- w /2 )]sin[^{y, + w/2) ], is 

electron mass, fc„ = {2mc/h?){E — En), k„ ~ 
^(2me/h^){En~E), En - {hy2m,)in^7ryw^), E is in- 
cident energy, is sum over all the evanescent modes. 

Explicit calculations of LDOS:We derive LDOS explic- 
itly from the internal wavefunctions for unit incident flux 
by using 



piE,x = 0,y = yi) ^ 



^ hVa 
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Where ipa{x,y) is the wavefunction when incident elec- 
tron is in the ath channel. It can be taken from Ref. Q 
and can be further simplified to give 
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In Eq. 8, sum over a runs over the propagating modes 
only while that for m is for all modes. 

LDOS from BTP formula: We consider two propagat- 
ing channels to present our results but we have checked 
that the results are the same for any number of channels 
including just one channel. For this case, Eq. 3 after 
summing over a and (3 gives 
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where primes mean derivatives wrt T^(0, yi). 



(10) 



V. VERIFICATION OF BTP FORMULA 

The derivation of Eq. 3 assumes that a small per- 
turbation to the actual system, allows us to expand a 
scattering matrix element as 

S%{E, V{x,,y,) T 5V{x,, y.)) = So.0{E, V{x,,y,)) 



j dx'M[5S^p{E, V{x[, y[))/5V{x[, y[)]5V{x[,y[) + . .(11) 



If we make a substitution of the type 

5 



5V{xr^y[) 



d 

dE 



(12) 



then we derive FSR. But this substitution being an ap- 
proximation, the FSR is also approximate. 

First of all, note that the role of evanescent modes on 
the derivative of the scattering matrix elements is that it 
renormalizes Sj according to the following relation. 



N 



(13) 



where gee = -^sin'^[^{yt + 'w/2)]. Note that the sum 
J2e 2k' ^® ^ converging series. It diverges as log[N] 
as iV ^ oo jS] . Here TV is the total number of evanescent 
modes. And so derivative w.r.t 7 will not exist for any 
arbitrary number of modes. Thus the expansion required 
to derive BTP formula is not defined at all energies and 
thus the concept of Larmor precession time fails. That 
means LDOS cannot be defined in terms of S matrix. 
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We shall also show that the expansion in Eq. 11 is valid 
at the Fane resonance and also explicit calculations of 
LDOS prove the correctness of BTP formula only at the 
Fano resonance. For this note that at the bound state 
2k^ ~ ~7' irrespective of the number of evanescent 
modes. Thus at the Fano resonance, convergence exists 
and the expansion in 11 holds good. Which means that 
at the Fano resonance, the BTP formula will give the 
correct delay time and hence also the LDOS. We verify 
this explicitly by numerically calculating the RHS of 8 
and 10. For example, let us truncate the series at the 
3rd term, that is we consider 2 evanescent modes only. 
In that case we show in Fig. 3 that the BTP formula 
is accurate. To calculate the derivative we have taken 
(57=0.001. However this kind of agreement does not occur 
for example when we consider 5000 evanescent modes 
with (57=0.001 (see Fig. 4). This can be cross checked 
by using Eq. 13. Note that now although the two curves 
do not coincide with each other exactly, they do coincide 
exactly at the energy corresponding to the bound state 
(i.e., where both curves peak). Of course one can take a 
smaller value of 5^ and get a better agreement between 
the two curves, but then again the same disagreement 
will be there if more evanescent modes are considered. 
Since the sum is not a converging sum, this will be a 
never ending story. One can see that = "^^"f ^ do 
not exist for an infinite number of evanescent modes as 
^ diverge as log[A^]. 

It is possible to provide further analytical arguments 
in support of our result. The RHS of Eq. 10 can be 
simplified to give 
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It was also shown in ref. l22l 
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This means that the functional derivative in 15 can be 
calculated by expanding the RHS of 16 and considering 
only the linear term in SV{0, jji), provided all higher order 
terms are finite. In the regime (1 + ^^) > V„ ^r^, 
this expansion can be done using (4) and (6), to give 
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It can be seen that the coefficients of higher order terms 
as well as the coefficients of linear terms in Sj, diverge as 



log[N] as well as higher powers of log[N], which implies 
that the functional derivative does not exist. 

Wang et al has pointed out some regimes where 
the BTP formula is violated due to the lack of gauge in- 
variance. They also calculated the correction terms in 
those regimes. However, they overlooked the violations 
that can arise due to non-existence of functional deriva- 
tives, and we point it out in this paper. Only way out is 
to truncate the number of evanescent modes to keep the 
error within acceptable limits. 



VI. GENERALIZATION 

The scattering matrix S of an extended system can be 
written as 5* = S*! (g) 5*2 (g) 5*3, where S2 is the scat- 
tering matrix of an infinitesimal region at r. Si and S'3 
are the scattering matrices of the regions to the left and 
right of r, respectively. Hence ^2 is the scattering matrix 
of a (5 function potential that we have to further perturb 
infinitesimally to see the changes in S2 ■ AH the change in 
S will be due to this change in 82- Thus complications of 
the BTP formula will depend on its complications for a 6 
function potential. If we want to integrate LDOS over an 
extended region then one has to take S function poten- 
tial like perturbations at many many places of the region 
and sum the changes they produce on the S matrix. The 
error will be added. 



VII. WDT AT FANO RESONANCE GIVES 
CORRECT SUPERLUMINALITY 

We have shown in section 5 that the BTP formula is 
exact at the Fano resonance. So the LHS in 12, operating 
on Saf3 will give the correct delay time as will be observed 
in an experiment. The substitution in 12, is exact at the 
Fano resonance follows from earlier results 0| that the 
FSR is exact at Fano resonance. Hence the RHS of 12, 
operating on Saf3 will also give the correct delay time 
as can be observed in an experiment. Hence, in spite of 
all non-stationary phase behavior, the WDT will exactly 
correspond to physical delay times at the Fano resonance 
and the negative slopes in Fig. 2 do mean that there is 
strong superluminality that can be observed. 



VIII. LARGE VIOLATIONS OF FSR DUE TO 
EVANESCENT MODES 



We see in Fig. 2 that away from the Fano resonance, 
there are large violations in FSR. Once we have explored 
the BTP formula we can now analyze the cause of it. 
First of all we would like to state that we are dealing with 
a system that is coupled to reservoirs (a grand canonical 
system) where the reservoirs can inject charge or absorb 
charge. So charge is not conserved, p — po can be ar- 
bitrary, and it may not be possible to relate it to 
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very strictly. But of course, as in all grand canonical sys- 
tems, at equilibrium charge and energy are as good as 
conserved and we can talk of Friedel sum rule. Eqn 12 
suggests that if there is a sample connected to semi-infite 
quasi ID leads then 

do 

= n[p{E) - po{E)]giobai = 



^ hVa 



mr)\^-\Mr)ndr 



(18) 



where r represents coordinate. Now the sample that ex- 
tends from — Ts to +rs is a grand canonical system and 
one can show that when the leads are sigle channel then 
(in the absence of evanescent modes) 
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i(E-Ei) error due to the substitution in Eq. 12. 

It depends on parameters like Ei that depends on the 
internal details of the potential. For multichannel leads, 
all the EnS, appear in the correction term. This amount 
is a negligible amount for practical purposes and the sub- 
stitution in 12 can still be used without involving large 
errors. We show below that in the presence of evanes- 
cent modes, there will be larger errors when we make the 
substittuion in 12. 

Ref. [23| did some explicit calculations to show 

[p{E)-po{E)\giobai ^'^-r — / propagating modes 

nVa J-oo 
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that gives 
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Here Va — hka/m^. denotes sum over all propa- 

gating modes and denotes sum over all evanescent 
modes. Taa = l^aale"''" is the ath diagonal element in 
the S'-matrix. t^^, is the transition amplitude from the 
ath propagating mode to the eth evanescent mode. 

So now if we follow the scheme of 19 then p{E) — 
Pq{E)—0 as 0. But if we consider 



[p{E)-po{E)] = J2 



hVa 



E 



(22) 



then we get the dotted curve in Fig 2. The first term on 
the RHS of 21, gives the charge lost to the reservoir. One 
can do the integration to find dxCos{2kaX + rja) = 
T:Cos{-qa)6{ka). So for ka > 0, charge is fully conserved. 
The difference in the dotted and dashed curves in Fig. 2, 
is the error due to the substitution 12. The error is still 
4(^eF-Ei) '-'^ single channel leads and this error is 

due to the substitution in 12. 

Substituion 12 means that change in scattering matrix 
element due to infinitesimal increase (decrease) in inci- 
dent energy, is the same as change in scattering matrix 
element due to constant infinitesimal decrease (increase) 
in potential over entire space. And also the constant de- 
crease in potential over entire space can be integrated as 
small changes locally and adding up the effect due to all 
such local changes. 

If we could have chosen a SV{r) at r < — rs that dis- 
turbs the propagating modes at r and does not disturb 
the evanescent modes at r then the error in substitution 
12 would have been as negligible as it is in ID. This is 
because although the evanescent modes exist in the leads 
they are to be considered as not existing in the leads but 
existing in the sample (this is clear from 19, 20, 21, 22). 
Indeed we see that in Eq 6, if we ignore any change in 
the second term on the RHS (which is due to evanescent 
modes) then the scattering matrix elements look exactly 
like that of a delta function potential in ID. But once we 
include the evanescent modes on the RHS of 6, we get 
diverging terms in the scattering matrix if we change E 
or V, and the substitution described in the above para 
involves much greater error as compared to ID. 



IX. SUMMARY AND CONCLUSIONS 

The BTP formalism is very crucial to understand 
mesoscopic transport beyond the Landauer conductance 
formula, that is beyond the linear and DC response. So 
far it has been verified for simple ID systems |2J| or sys- 
tems where there are localized states completely decou- 
pled from the leads !25| . Violations of BTP formula can 
also occur in situations where there is absence of gauge 
invariance |2^. We have shown that BTP formalism is 
not exact in Quasi-one-dimension (QID) due to the pres- 
ence of evanescent modes. These evanescent modes make 
the scattering matrix singular in nature and the series 
expansions in Eq. 11. required to show the equality 
between p and breaks down. If there are only a few 
evanescent modes in the system, then the formula may 
be acceptable for practical purposes. For extended poten- 
tials and large number of evanescent modes, it may not 
be practical. It definitely cannot be used as a definition 
for DOS. DOS has to be defined in terms of the internal 
wave function and Hamiltonian. Scattering matrix will 
not contain all the information. We have also proved 
that in spite of all non-stationary phase effects, WDT 
{-^arg(taa)) correctly give the superluminal times and 
there is superluminality at the energy where the phase 
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drops of transmission amplitudes occur (for example at 
Ea—85 in Fig. 2). Experimentalists have always tried 
to find situations wherein they can create larger super- 
luminality. A delta function potential in a QID wire 
can create limitless superluminality. So far, experiments 
to observe superluminality, only considered systems with 
Breit-Wigner type resonances pq|. 

Finally in this paragraph, we also make some com- 
ments on possible future research, that are not related to 
the theme of this work. Features of superluminality in 
case of Fano type resonances should be explored exper- 
imentally. Also as it seems that semi-classical formulas 
like WDT and FSR are valid at Fano resonances, in spite 
of strongly energy dependent scattering phase shifts, it is 
possible that an incident wave packet is not dispersed by 



scattering at Fano resonance. Parts of it will be trans- 
mitted to different channels without shape distortions. 
Rather, if wavepackets are simultaneously incident on 
the scatterer from different channels, then after scatter- 
ing the outgoing wavepackets in different channels will 
be similar to the incoming wavepackets. Their centroids 
may remain unchanged before and after scattering. This 
normally happens for solitons and very rare in quantum 
mechanics. 
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Figure captions 

Fig. 1. Here we show a quantum wire of width w 
with a delta function potential V{x,y) = ^5(x)5(y — yi) 
situated at x. We consider scattering effects when the 
incident electron is from the left. The sub-bands on the 
left of the impurity is denoted as 1 for the first mode and 
2 for the second mode. Similarly the sub-bands on the 
right of the impurity is denoted as 3 for the first mode 
and 4 for the second mode. 

Fig. 2. The dotted curve gives G=n{p{E) - po{E)), 
where p{E) is DOS at energy E in presence of the scat- 
terer and po{E) is DOS at energy E in absence of the 
sctterer. The system considered is in Fig. 1, and G 
is in units of g = (mew'^/h?). The dashed curve gives 
Q — d9f/dE in units oi q — g. The solid curve gives 
arg{tii) in radians, displaced on the y-axis by 0.1. The 
dash-dotted curve gives arg{t22) in radians, displaced on 
the y-axis by 0.8. We have taken 7 = —15, 7/^=0.45 and 2 
propagating modes along with 2 evanescent modes. The 
X-axis is energy in units of a = w'^. 

Fig. 3. The system under consideration is shown in 
Fig. 1. with two propagating modes. Solid line gives 
H — p{E,0,yi), that LDOS calculated from internal 
wavefunction, in units of h — {2mw/h^). Dashed line 
gives J — i/{E,0,yi), i.e., LDOS calculated using BTP 
formula, in units of j = h. All parameters are the same 
as that considered in Fig. 2. The x-axis is energy in units 
of a — w^. 

Fig. 4. Here we consider 5000 evanescent modes and 
7 — —1.5. Everything else is the same as in Fig. 3. 
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